Topological Wilson-loop area law manifested using a superposition of loops 
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We introduce a new topological effect involving interference of two meson loops, manifesting a 
path-independent topological area dependence. The effect also draws a connection between quark 
confinement, Wilson Loops and topological effects. Although this is only a gedanken experimment 
in the context of particle physics, such an experiment may be realized and used as a tool to test 
confinement effects and phase transitions in quantum simulation of dynamic gauge theories. 



I. INTRODUCTION 

Topological and geometric effects are fundamental 
quantum-mechanical phenomena. They appear in var- 
ious physical contexts, as the Aharonov-Bohm and 
Aharonov-Casher effects [TJ [2], Berry's phase [3] and 
other models. Such effects are manifested by accumu- 
lated topological or geometrical path-dependent phases, 
which are observed in interference experiments. 

In quantum field theory and particle physics, a topo- 
logical phase similar to the AB effect appears in the con- 
text of the Wilson loop operator [H [5] , which is an or- 
der parameter manifesting the disorder in the confining 
phase of a gauge theory. The Wilson loop operator along 
a curved spacetime path C is 

W(C) = p( e ^c A ^\ (i) 

where P stands for path ordering (see [B], for example). 

Wilson Loops manifest confinement through the area 
dependence of their expectation value, and thus are an 
important test for confinement. However, they involve 
a product of operators along a loop and thus are non- 
local. Furthermore, they can be interpreted as transition 
amplitudes, and thus their phases do not contribute to 
the related probabilities. 

In this paper, we present a method to observe the area- 
law manifested in confining theories using superposition 
and interference of mesons, which unlike in the Wilson 
loop approach, contains the relevant phase as a relative 
one. We draw the connection between our method and 
the well-known Wilson Loop. 

The paper is organized as follows: first, in section II, 
we consider the interference effect of two mesons, consist- 
ing of static quarks, in a superposition, and show how to 
gain the string tension from it. The relation to the Wil- 
son Loop operator is drawn. In section III, we allow one 
of the quarks to be dynamic, modeling it as a particle 
in a moving harmonic potential, and show how to obtain 
the string tension in that case, using an exact solution 
of Schrodinger equation. Finally, in section IV, we dis- 
cuss the possibility of realizing the proposed gedanken 
experiment using a quantum simulator [7]. 



II. SUPERPOSITION OF LOOPS: STATIC 
QUARKS 

Free quarks can not be found in nature, but rather form 
hadrons, due to the phenomenon of quark confinement 
[3]. A quark and an anti-quark, attached to each other by 
a confining flux-tube, form a meson, which is the simplest 
hadron in QCD. The static potential between the quarks, 
as a function of their distance R, takes the form 

V (R) = jR (2) 

for large values of R, where 7 = 7(g 2 ) is called the string 
tension, and g is the coupling constant pll IMTU]. 

Consider two static quarks, initially separated by dis- 
tance xq. By static we mean that we treat them as exter- 
nal sources of electric field, which do not have their own 
dynamics, but can be externally moved. Taking con- 
finement into account, we know that the two quarks are 
connected by a long flux-tube, with length R, forming 
a meson (Quarkonium). If we allow one quark to move 
along some trajectory x (t) — xo+d (t), we get the Hamil- 
tonian 

Hstat = 7 (*o +<*(*)) (3) 

where we generally assume that 7 is an unknown quan- 
tity. 

Next, assume that the moving quark has some two 
internal energy levels s, which we denote by |t) , \i) - 
eigenstates of the Pauli Matrix a z . We define the pro- 
jection operators to the subspaces of internal levels as 
P-f = |t) (t| 7 Pi = \i) (4-1, an d introduce an internal level 
dependent paths d s (t), i.e., 

H s tat = 7 (x + J2 d s (t) Pa) (4) 

We wish to consider the interference effect of two 
mesons in superposition, varying their length. We shall 
consider a superposition of a meson with a f fermion and 
a meson with a 4- fermion, initially with the same length 
(i.e., d s (t) = for t < 0), which are stretched by moving 
the right fermion in two opposite directions, and then, 
at time T > 0, brought back to the same length (i.e., 
d s if) = also for t >T. The opposite directions impose 
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FIG. 1. The gedanken-experiment. Two fermions qi,q2 form 
a meson of length xo at time t = 0. While qi remains static, 
§2 is externally moved (with no dynamics of its own) in a su- 
perposition of two possible trajectories, corresponding to its 
internal states |t) and creating a superposition of mesons. 
The path of the f fermion is drawn in red, and the path of 
the J, one - in blue. At time t = T, both the trajectories ar- 
rive again at xq, and the state is then mixed using a unitary 
transformation (equation |7|), causing an interference effect. 
In case of confinement, the interference phase depends on the 
area A (defined in equation (JsJ) ) in case of confinement - as 
does the Wilson-Loop. The phase is accessible from the prob- 
abilities to find q2 in either of its internal states (equation 




more conditions on the paths: (t) > 0, (t) < 0, and 
Since the .J, fermion should not "go through" the static 
fermion at x = 0, we also demand |d^(t)| < x . The 
initial state is, of course, 

hKo)) = |t*> = ^(lt) + |;» (5) 

Solving Schrodinger equation, one gets that 

\if,(T)) = (er 11 ^ <M* ) dt \\) +e -^S a d v (t )dt |^ 

V2 V 

(6) 

and so in order to measure 7, one has merely to cause 
an interference between the two states. This is done by 
applying the rotation U = e -1 !" 5 on the state at t = T, 
and one obtains, up to a global phase, 

UH{T)) = (sin (?f) |t) +*cos (?f) II)) (7) 

where 

A= f {d^t^-d^t'^dt 1 (8) 
Jo 

is the area enclosed between the two trajectories (as in 
figure [lj. Thus, the probabilities to find the system on 
each of the internal levels are 

P t = S i^(^);n=co^(^) (9) 



By performing such an interference experiment, and 
measuring the phase, knowing the area difference one can 
calculate the string tension 7. Moreover, if the phase 
does not exhibit such an area law, it means that the 
system is not within a confining phase, and hence this 
measurement can be used for probing the phase as well. 

III. DYNAMIC QUARK MODEL 

Next we wish to introduce a simple model of dynam- 
ical charges, in order to examine the corrections to the 
latter static case. In particular, we wish to examine de- 
coherence and destructive interference due to excitations 
of the mesons. For the sake of simplicity, we assume that 
one of the quarks is static (the one placed in x = 0), 
and that the other one, having a mass m, is trapped in 
a harmonic potential, centered around x = xq > 0. The 
Hamiltonian of the system takes the form 

H = - h-w 2 (i-i ) +jx (10) 

2m I 

Assume that the dynamic quark has two internal en- 
ergy levels, denoted and treated as before. We introduce 
an interaction between the internal and external degrees 
of freedom, of the form 

H int = xY^G s {t)P s (11) 

s 

Where G s (t) are opening functions which are zero for 
times t < 0, t > T, and are assumed to be smooth 
enough, i.e. at least their first and second time deriva- 
tives vanish for t = 0, T. 

Define d s (t) — — m 3 ^ 2 ; (the same conditions on d s (t) 
apply as in the static case, of course, and that poses con- 
ditions on G s (t)), x s (t) — xq — + d s (t). Then one 
gets, after completing the square, that the total Hamil- 
tonian is (neglecting constants) 

H = ^2(H,+fd.(t) + f B (t))P a (12) 

s 

where the first part is just a Harmonic oscillator, with the 
center of its potential following classically the trajectory 
x s (t): 

H s = ^ + \rnuj 2 (x - x s (t) f (13) 
2m 2 

and/ s (i) = G s (t)x Q -^Gl (<). 

A. Exactl Solution of Schrodinger Equation 

The next step is the solution of Schrodinger Equation 
for our Hamiltonian: 

l ^U)=H\^) (14) 
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Noting that the internal levels are not changed by the 
Hamiltonian, we can first solve for a given internal level 
s and consider only the harmonic part. Thus if we set 
for a given s, 



^(»)\ = e -<(7/ '«».(*')*'+/o /.(*')*') $w 



(15) 



t/>( s ) ) is the solution of a Schrodinger equa- 



we get that 

tion for a classically moving Harmonic potential: 
. d 



dt 



(16) 



using the solution of this equation in x-space, we get that 
the solutions (not eigenstates) are 



^ (X,t) = e ^ 8) (*)e-^*e-^/o d.(t>)tfjm(d.(t)+q.(t))( X -x.W) Xn {x _ ^ {f) _ qg (t)) 



(17) 



where E n , \ n (x) are the energies and eigenstates of a 
"regular", fixed-potential harmonic oscillator, (t) — 
- /„' f s (f) dt' + f /„' ell it) f J* (<g (t) - u\l (t)) and 
<j' s + cj 2 ^ = -d s . 

Our initial condition is (x,0) = xa (x — x s (0)) = 
Xo {x — xq). For t < 0, the oscillator is supposed to be 



in its non moving ground state, and thus we expect that 
q s (t) = for these times. Thus q s (0) = 0. Using the 
continuity of the e.o.m. of q s , we get that q s (0) = as 
well. From the smoothness of the opening functions we 
know that d s (0) = 0. Thus we conclude that in order to 
start from a local ground state, the solution must be 



tP^ (x,t) = e ^ s) (t) e -^ot e -n/o t ' i =(*') d *' e lm (^( t )+^( t ))( a: - a; =(*)) Xo ( x - Xs (t) - q s (t)) 



(18) 



Next, we wish to interpret this solution in terms of lo- 
cal instantenous eigenvalues. That is, the states \n(t)), 
defined as eigenvalues of H s (t) in the "frozen" time t. 
We already know that the system starts at t = with 
the eigenstate |0(0)), but what's later? In order to do 
that, we define y = x — x s . Writing H s (t)in terms of 
y at a fixed t yields this diagonalization: it is merely a 
transformation to a frame which moves with the poten- 
tial, which is its instantaneous rest frame (IRF). Consider 
the position-dependent part of (x,t), in terms of the 
IRF. Call it <j>^ (y,t): 



(t) 



im{d s (t)+q B (t))y -iq s (t)p y 



0(t)> (19) 



in the IRF basis, this is a coherent state - Poissonian 
distribution of \n {£)) states: 



where 



a s (t) 



(t)} = e i"l(d s (t)+q B (t))q s (t) ^ ( 2 Q) 

n (d s (t) + q s (t)) 



moj 



q s (t) 



V2 



(21) 



B. Superposition and Coherence 



Suppose we start, at t = 0, with an initial state 

|^(0)) = |0(0)>|tx> = ^|0(0))(|t) + |4» (22) 



4 S) (*) 



then, defining 
\m [d s (t) + q s (t)j q s (t) and using the solutions 
from the previous section, we get 



WW) 



ot/2 



(23) 

Let us understand the meaning of this state. We start 
at t = with a superposition of two states with two 
different values of s. Both of them are in the ground 
state of an oscillator, centered around the same position 
(x + (0) = x_ (0) = xo). Then we "move" the wavefunc- 
tions together with the potential: the opening functions 
are translated to the trajectories d s (t). Each element 
of the superposition "goes" through another path, since 
the moving of the potential depends on the internal level. 
A superposition of coherent states in the terms of the 
local IRFs is created, and eventually, at t = T, both 
the interaction functions are closed, i.e. the two poten- 
tials experienced by internal levels overlap again, and 
x+ (T) = x- (T) = x . 

Next, trace out the oscillator degrees of freedom at 
t = T, to obtain an internal-level density matrix. The 
density matrix at t = T is 



P = 



M 1 

2 1 e -,(* m -* ( «) B * 



B 



(24) 



where $W = % s) (T) - 7 J* d s (t') dt' and B 
e -Kl a t(T)| 2 +K(T)| 2 ) a t (TK(T) 
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One can see that the phases depend on two parts: one, 
CO' i s totally calculable. The other is 7 dependent 
and can't be calculated unless 7 is known. The 7 de- 
pendent phases are only global (the probabilities - the 
diagonal terms in p do not depend on them). In order 
to observe the phases, we wish to cause an interference, 
i.e. to rotate the state: We act on the system at t = T 
with the rotation operator U = e~ l ^ a ' v . One can use the 
diagonal terms in the new density matrix UpW to de- 
termine the probabilities to measure each of the internal 
levels. The probabilities are 

P t = i (1 - cos (A$) RcB + sin (A$) ImB) 

Pj, = ~ (1 + cos (A$) RcB - sin (A$) ImB) (25) 

where A<i> = $(t) — $(40 We see that now the probabili- 
ties depend on the phase difference, and that it became a 
relative phase indeed. However, how do we use it in order 
to measure the string tension? First, one must note that 
B is governs the visibility of the interference: if B = 0, 
one gets equal probabilities to both measurement out- 
comes, and hence no information can be gained and the 
interference is lost. We shall consider the effect of B in 
detail. However, let us first focus on the role of the phase 
difference, assuming it is not screened by B. 

The phase difference consists of two parts. The first 
one, $ = $o T) (T) - $^ (T), is 7 independent, and 
once the phase is obtained from the probabilities, it can 
be subtracted. Hence we are left with 



7/ {<h(t')-d x (t'))dt' = - 7 A 



(26) 



where A is the area enclosed between the two paths 
in spacetime! Exactly the area dependence which is 
expected in any abelian and non-abelian gauge theory 
within the confining phase.. Thus, the interference ef- 
fect measures the string tension 7 in case of confinement; 
Otherwise, the area law will not be manifested and thus 
being outside the confining phase can be probed this way 
as well. 

Let us discuss the effect of B. In order to understand it, 
we calculate it explicitly; Using the definitions of B,a s , 
we get that B = B^e 15 , where 



B Q = e 



= p-K I T(9t(T)+94(T)) 2 + S(<?tm+'34m) 2 ) 



S=-(q t (T) qi (T)-qUT)q t (T)) (27) 



The interference is maximal when 5 = 1; That cor- 
responds to a = 0, which means that the final state is 
an eigenstate rather than a Poissonian superposition. In 
that case, q s (T) = q s (T) = 0. Then we get 



D (B=1) • 2l^0~lA\ D (B=1) o 

PI = sin ; P, — cos 



$0 -7^ 
2 

(28) 



If we wish to consider the cases in which there is some 
disturbance to the interference, but it is negligible, we 
should consider the limit B — > 1. This is obtained when 
the conditions yf^fq s (T) < 1, y/^q a (T) < 1 are met. 
This can be understood in terms of uncertainty principle: 
the final coherent state is displaced in phase space. The 
displacement in x is q s (T) , and the first condition is met 
if we require it to be much smaller than the ground state's 
Aa;; The displacement in p is mq s (T), and the second 
condition is met if we require it to be much smaller than 
the ground state's Ap. 

Let us see what is the limitation on the trajectories 
d s (T) , if one wishes to get a good interference according 
to this criterion. In order to do that, let us write the ex- 
plicit solution for q s (t). We wish to solve the differential 
equation q s + uj 2 q s = —d s . Previously, we have obtained 
the initial conditions q s (0) = 0, q s (0) = 0, and thus the 
homogenous solution is zero, and we are left only with 
the particular solution, which can be found using Green's 
function: 



q s (t) = / sin (u (t - if)) d s (t 1 ) dt' 



q s (t) = - cos (ui (t - if)) d s (*') dt' (29) 
Jo 



Demanding \q s (T)\ -C Aa; = 7==, we get a condition 



of the maximal accelaration a„ 



(<*-(<)): 



< 



1 

T V 2m 



(30) 



and the very same condition is obtained from demanding 
m\q s (T)\ « Ap = y/^f. Thus we conclude that the 
interference is not ruined if the charges' accelerations are 
small enough all along the paths. 

Finally let us examine the relation to the static model 
of the previous section. If one wishes to consider the case 
of two static quarks, since the initially dynamic quark is 
not free but rather subject to a harmonic potential, tak- 
ing the m — > 00 limit wouldn't be enough, and u> — > 00 
should be considered as well, keeping their ratio constant. 
That will result in a highly localized quark, with a very 
large momentum uncertainty. However we do not mea- 
sure the momentum and thus there is no problem with 
that, the constant ratio will assure us that we can still 
use a nonzero maximal acceleration satisfying condition 
(301, and the interference will not be lost. There is a 
slight difference resulting from <I>o, but it can be sub- 
tracted in the results, or taken into account by adding 
some time-dependent terms to cither of the Hamiltoni- 
ans (static or dynamic) - then a full equivalence between 
the model can be reached in the static limit. 
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IV. RELATION TO THE WILSON LOOP 
OPERATOR 

Next, let us examine the relation of the proposed 
method with the Wilson Loop operator approach. 

The area A , enclosed between the two paths of the 
two quarks in superposition, is the same area on which 
the Wilson-Loop of a single quark, moved in spacetime 
along the union of the paths, would depend, and with 
the same string tension. Thus, the interference phase 
and this Wilson-Loop's phase are the same. We shall 
now describe how to derive a quantitative relation to the 
Wilson Loop operator in the case of static quarks and 
strong coupling limit. 

Denote the state of heavy (static) QQ, separated by 
distance R by \R (g 2 )), and write it in terms of the gauge 
field degrees of freedom. Assuming confinement, this cor- 
responds to a quarkonium state, where the two quarks are 
connected by a flux tube, and thus 

i/|i?(.g 2 ))= 7 (.9 2 )i?|i?(.g 2 )) (31) 

where H is the Hamiltonian of the gauge field; In partic- 



ular, in the strong coupling limit one gets 

lim \R (g 2 )) = P (e l ^ A ^\ \vac) (32) 

where P stands for path ordering, as, for example, in [6]). 

We wish to calculate, within the strong coupling 
limit, the expectation value (in Minkowski space) of the 
Wilson-Loop operator, corresponding to the loop de- 
picted in figure (|]a). This can be decomposed to four 
different parts: I, II, III and IV. The IV part contri- 
bution is zero if we work in the temporal gauge. Let us 
now see the contribution of the other three parts of the 
loop. In order to do that, let us discretize the function 
d (t) as in figure ^p) . Define 

F(x l ,x 2 ) = p(e if *? A » dx ") (33) 

acting on the vacuum, this operator creates (in the strong 
limit) a flux tube from x\ to xi\ i.e., the static charge 
in x\ is raised and the static charge in X2 is lowered. 
Using these terms, the Wilson Loop operator we wish to 
calculate becomes 



l 

W (C) = e lHT F (x , 0) e~ lHT ]J (e im «F {x n . x ,x n ) e~ iHt -) F (0, x )) (34) 

n=N 

since t n+ i — t n = £ = AT, this expression simplifies to 

l 

W (C) = e iHT F 0r o , 0) J] (e~* HAT F (x n ^,x n )) e~ iHAT F (0, x )) (35) 



Finally, let us calculate the expectation of the Wilson 
Loop in the vacuum state. Taking the vacuum energy 
as zero, the left e lHT contributes 1. The right F (0, Xo), 
acting on the vacuum, creates a flux tube between x — 
and x — xO, and the e~ tHAT on the left of it contributes 
a phase of e" !7I ° Ar since it is an eigenstate. Then, each 
of the F (x n -i, x n )'s in the product shortens or stretches 
the flux tube, and the e ~ tHAT to the left of it adds up 
a phase of e~ 1/yXn , This is a process of creating a 
flux tube with length Xq at t — 0, changing its length 
according to d(t), until t = T where its original length 
is regained and it is destroyed. The amplitude for this 
process is thus, according to the given explanation, the 
expectation value of the Wilson Loop, and it is 

N-l 

(W (C)) = e "=» (36) 

Taking back the continuum limit, we take N — > oo, or 
AT — > 0. This transforms the sum to an integral, and 
this integral is equal to A, the area enclosed by the curve 



- the area law of confinement: 

(W (C)) = e-^fi *(*')<"' = e -^ A (37) 

Thus we can conclude, that in the strong coupling 
limit, taking the solution from equation (|6|, 

(s\iP(T)) = -±={W(C s )) (38) 

(where C s denotes the curve enclosed by the motion of 
the sth element of the superposition) - so we see that 
indeed, in this limit, the transition amplitude is the cor- 
responding Wilson Loop, and it is a global phase, so in 
order to obtain knowledge about it one has to transform 
it to a relative one, using the above interference prescrip- 
tion. 



V. QUANTUM SIMULATION 

Quantum simulations [7j are a rapidly growing physical 
field, based on the idea that quantum systems can simu- 



() 




AT = T / N 



FIG. 2. (a) The phase of the Wilson Loop operator depends 
on the area enclosed within it, in confinement phase, (b) 
Discretization of the curve, as explained in the text. 



late each other. Thus, quantum systems which are inac- 
cessible for measurement, can be simulated using other 
quantum systems, which are controllable, accessible and 
measurable in the laboratory, such as cold atoms in opti- 
cal lattices [HJ[T2], trapped ions [T3j and other systems. 
These systems serve as an " analog quantum simulators" . 

Although this is a gedanken experiment, it can be re- 
alized in the context of quantum simulations of confining 
phases of dynamical gauge theories [21 [TS] . In such sys- 
tems, one could use lasers to create a superposition of 
fermions, which results in a superposition of mesons, and 
to enable the mixing and interference required for the 
measurement |16) . 



interference of mesons. Interestingly, the two states in 
superposition experience different electric fields. This is 
a reminiscent of the concept of "private potential" |17j . 
By exploiting an area-dependent phase due to the linear- 
ity of the static quark potential in the confining phase 
of gauge theories, we have observed confinement using 
an interference of the mesons. This allows to convert 
the global phase (which appears, for example, in Wilson 
Loops which are related to transition amplitudes), to a 
relative phase, observable in probabilities. 

In the Coulomb phase (as in 3+1 QED) or in any other 
V oc RP phase, with (3 ^ 1, the appropriate gauge field 
state does not include a flux tube as in the confining 
phase. Hence, the final probabilities will not manifest a 
simple area dependence. This can be used to probe a 
transition between confining and non-confining phases. 

Although it is only a gedanken experiment in the con- 
text of particle physics, such an experiment, in its lattice 
version, may be realized and used to observe confinement 
effects and phase transitions within a quantum simula- 
tion of confining gauge theories [21 [TS] . 
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